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We introduce the concept of bond spin current, which describes the spin transport between two 
sites of the lattice model of a multiterminal spin-orbit (SO) coupled semiconductor nanostructure, 
and express it in terms of the spin-dependent nonequilibrium (Landauer-Keldysh) Green functions 
of the device. This formalism is applied to obtain the spatial distribution of microscopic spin 
currents in clean phase- coherent two-dimensional electron gas with the Rashba-type of SO coupling 
attached to four external leads. Together with the corresponding profiles of the stationary spin 
density, such visualization of the phase- coherent spin flow allow us to resolve several key issues for 
the understanding of mechanisms which generate pure spin Hall currents in the transverse leads 
of ballistic devices due to the flow of unpolarized charge current through their longitudinal leads: 
(i) while bond spin currents are non-zero locally within the SO coupled region even in equilibrium 
(when all leads are at the same potential), the total spin current obtained by summing the bond 
spin currents over an arbitrary cross section is zero so that no spin can actually be transported by 
such equilibrium currents; (ii) when device is brought into nonequilibrium steady current state by 
external voltage difference applied between the longitudinal leads, only the wave functions (or Green 
functions) around the Fermi energy contribute to the total spin current through a given transverse 
cross section; (iii) the total spin Hall current is not conserved within the SO coupled region; however, 
it becomes conserved and physically well-defined quantity in the ideal leads where it is, furthermore, 
equal to the spin current obtained within the multiprobe Landauer-Biittiker scattering formalism. 
The local spin current profiles crucially depend on whether the sample is smaller or greater than 
the spin precession length, thereby demonstrating its essential role as the characteristic mesoscale 
for the spin Hall effect in ballistic multiterminal semiconductor nanostructures. Although static 
spin-independent disorder reduces the magnitude of the total spin current in the leads, the bond 
spin currents continue to flow through the whole diffusive 2DEG sample, without being localized as 
edge spin currents around any of its boundaries. 

PACS numbers: 72.10.Bg,73.23.-b,72.25.Dc,71.70.Ej 



Recent experimental observation of the spin Hall ef- 
fect 1,2 opens new avenues for the understanding of fun- 
damental role which spin-orbit (SO) couplings^ can play 
in transport and equilibrium properties of semiconductor 
structures. While SO coupling effects are tiny relativistic 
corrections for particles moving through electric fields in 
vacuum, they can be enhanced in solids by several orders 
of magnitude due to the interplay of crystal symmetry 
and strong crystalline potential. 3,4 Furthermore, harness- 
ing of spin currents induced by the spin Hall effect offers 
new possibilities for the envisioned all-electrical manipu- 
lation of spin for semiconductor spintronics applications 5 
where electrical fields can access individual spins on short 
length and time scales. 

The principal macroscopic phenomenological manifes- 
tation of the spin Hall effect is unique: the transverse 
spin current, which is pure in the sense of not being 
accompanied by any net charge transport in the trans- 
verse direction, emerges as a response to conventional 
unpolarized charge current in the longitudinal direction 
through a paramagnetic system in the absence of any ex- 
ternal magnetic field. When such current hits the sample 
boundary, it will deposit nonequilibrium spin accumula- 
tion&L2i2iiS at the lateral edges of the sample attached 
to two longitudinal electrodes, 9,10 as detected optically 
in recent breakthrough spin Hall experiments. 1,2 

However, there are several apparently disconnected 



mechanisms capable of inducing the spin Hall currents. 
Nonetheless, they share the necessity for some type of SO 
interaction which couples the spin and charge transport. 
For example, impurities with SO interaction will deflect 
spin-j (spin- j) conduction electrons predominantly to the 
left (right) in the scattering process, thereby generating 
the extrinsic transverse spin Hall current. The theory 
of the extrinsic spin Hall effect has been around for sev- 
eral decades and it has recently been revisited 11 to 
argue its major role in one of the two recent seminal ex- 
perimental observations A However, the extrinsic effect, 
which crucially relies on the presence of impurities with 
skew-scattering and does not involve any SO coupling in- 
duced modification of the quasiparticle energy spectrum, 
is a rather small effect (unless one invokes particular en- 
hancement mechanisms involving intrinsic SO coupling 
in the bulk crystal which contributes a SO term to the 
impurity potential 11 ) whose precise magnitude has been 
tantalizingly hard to estimate. 12 

Thus, a strong impetus for the revival of interest in 
the realm of the spin Hall effect (s) has ascended from 
recent predictions for substantial pure spin current flow- 
ing through: (a) infinite homogeneous SO coupled semi- 
conductors in the clean limit ^iiii where the strong SO 
coupling induces the spin-splitting of the quasiparticle 
energies; or (b) the electrodes of multiprobe finite-size 
mesoscopic nanostructuresi 5 ^^ made of such materi- 
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als. However, the theory of the intrinsic spin Hall effect 18 
in the bulk of infinite semiconductors is formulated in 
terms of the spin current density which is not conserved^ 
in a medium with SO coupling and, moreover, can be 
non-zero even in thermodynamic equilibrium. 19 Reexam- 
ination of such unusual features of the intrinsic spin Hall 
current have lead to two major controversies: (i) its de- 
pendence solely on the whole SO coupled Fermi sea and 
the spin-split band structure (i.e., the anomalous velocity 
due to the Berry curvature of the Bloch state s 13 ! 20 ) , with- 
out any connection to nonequilibrium distribution func- 
tion that characterizes standard charge transport in the 
longitudinal direction, 21 have prompted arguments that 
such currents do not really transport spin or induce spin 
accumulation that would be useful for spintronics appli- 
cations^!^ (ii) for linear in momentum SO couplings, 
such as the Rashba or linear Dresselhaus coupling 4 in a 
two-dimensional electron gas (2DEG), numerous pertur- 
bative analytical2iS4 and non-perturbative numerical ex- 
act diagonalization 25 studies find that bulk intrinsic spin 
Hall current (averaged over an infinite system^ ) vanishes 
for arbitrary small disorder, while being able to survive as 
edge spin current near the sample-electrode interfaces^ 
or for some intermediate values of the spin-independent 
disorder strength 28 (for which the sharp boundary be- 
tween the intrinsic and the extrinsic mechanisms becomes 
blurred2&) . 

On the other hand, the mesoscopic spin Hall current, 
predicted to flow out of ballistic phase-coherent samples 
of the SO coupled semiconductorsiSiASiAiiSliSiiSS through 
the attached transverse ideal (i.e., free of spin and charge 
interactions) leads, is conserved throughout the leads, de- 
pends only on the wave functions (or Green functions) at 
the Fermi surface (at low temperatures T — > 0), and it 
is resilient to rather large static disorder within the dif- 
fusive metallic regime^ii However, the theory of the 
mesoscopic spin Hall effect is formulated in terms of 
the Landauer-Biittiker-type transmission formalism 3 ^ for 
spin currents ^5iASi2£ which connects asymptotic scatter- 
ing states in the leads without requiring any informa- 
tion about the quantum-mechanical probabilities for spin 
and charge propagation between the two points inside 
the sample. Technically, to obtain the spin Hall current 
flowing through the leads of a mult iter minal device, one 
only needs the spin-dependent retarded real-space Green 
function connecting the sites residing in different leads 
and no information about its values between the points 
within the sample. 15 

Thus, many recent debates on the very existence of 
the spin Hall effect in ballistic SO coupled semiconduc- 
tor systems could be closed by visualizing the spatial 
details of the spin flow through experimentally accessi- 
ble Hall bridges — from the SO coupled sample toward 
to attached electrodes with no SO interactions. The 
analogous studies of the spatial distribution of charge 
flow were essential in understanding the nature of quan- 
tum Hall transport (bulk vs. edge 35 ) in mesoscopic Hall 
bridges Furthermore, recent advances in multifari- 



ous scanning probe experimental techniques have made 
it possible to go beyond conventional transport measure- 
ments of macroscopically averaged quantities and image 
phase-coherent charge flow through a single 2DEG in 
quantum Hall or quantum point contact regime where 
the host semiconductor heterostructures is subjected to 
high or zero external magnetic field, respectively^ 

In particular, this type of insight can be obtained effi- 
ciently within the framework of lattice models of meso- 
scopic devices and the corresponding bond charge cur- 
rents 34,39 which yield a detailed picture of the charge 
propagation between two arbitrary sites of the lat- 
tice^fiiii Here we provide a tool that makes it possible 
to obtain the spatial details of the spin flow on the scale 
of few nanometers by introducing the bond spin currents, 
which represent the analog of bond charge currents (and 
a lattice version of the spin current density). 

As shown in Sec. [IJ the bond spin currents can 
be computed efficiently in terms of the spin-resolved 
Keldysh Green functions 42 applied to the Landauer set- 
up where a finite-size SO coupled semiconductor sam- 
ple is attached to many semi-infinite ideal leads 
They are evaluated in Sec. |H] for a paradigmatic meso- 
scopic spin Hall generator — a four-terminal ballistic two- 
dimensional electron gas (2DEG) with the Rashba type 
of SO coupling — to show the spatial profiles of spin cur- 
rents and spin densities, thereby, revealing the features of 
the spin Hall transport on nanoscales. In Sec. lIIII we show 
that spatial distribution of the local spin current densi- 
ties remains non-zero throughout the whole 2DEG even 
in the diffusive transport regime, in contrast to some pre- 
vious conjectures2fi*2i where only the edge spin currents 
would survive the disorder effects in Rashba spin-split 
systems. The integration of the linear response bond 
spin currents over the transverse cross sections allows us 
to connect in Sec. lIVI the spin transport within the sample 
to the total spin Hall currents, which are obtained from 
the Landauer-Biittiker multiprobe spin current formulas, 
flowing in an out of the 2DEG through the leads as a re- 
sponse to the applied voltages at the device boundaries. 
We conclude in Sec. 



I. BOND SPIN CURRENTS IN 
MULTITERMINAL NANOSTRUCTURES: 
LANDAUER-KELDYSH APPROACH 

The conservation of charge implies the continuity equa- 
tion in quantum mechanics for the charge density p = 
e|\I/(r)| 2 associated with a given wave function \£(r) 



dp 

at 



Vj = 0, 



(i) 



from which one can extract the expression for the charge 
current density 



j = eRe [^ f (r)v^(r)]. 



(2) 
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This can be viewed as the quantum-mechanical expecta- 
tion value [in the state ^(r)] of the charge current density 
operator 



n(r)v + vn(r) 



(3) 



which follows from the classical charge current density 
j = en(r)v via quantization procedure where the particle 
density n(r) and the velocity v are replaced by the corre- 
sponding operators and symmetrized to ensure that j is a 
Hermitian operator^ In SO coupled systems j acquires 
extra terms since the velocity operator ifw = [r, H] is 
modified by the presence of SO terms in the Hamiltonian 
H. For example, for the effective mass Rashba Hamil- 
tonian of a finite-size 2DEG structure (in the xy-plane) 



H = + - (py& x - p x a y ) + V coni (x, y), (4) 



2m* ' h 
the velocity operator is 



P OL A 



(5) 



where e x and e y are the unit vectors along the x and 
the 2/- axis, respectively. Here p = (p x ,P y ) is the mo- 
mentum operator in 2D space, & = (a x , a y , a z ) is the 
vector of the Pauli spin matrices, a is the strength of 
the Rashba SO coupling 3,4 arising due to the structure 
inversion asymmetry, 4 - and V con {{x 1 y) is the transverse 
confining potential. 

In contrast to the charge continuity equation Eq. (QJ, 
the analogous continuity equation for the spin density 



A h 



[*t( r )ai*(r) 



dt 



V-J l = Sl, 



(6) 



(7) 



contains the spin current density 

^ = 2 W 2 (r) ' 

as well as a non-zero spin source 



Sl = ^Re (V(r)^[ff,a;]*(r)). (8) 

The non-zero S l s ^ term reflects non-conservation of 
spin in the presence of SO couplings. Thus, the plausible 
Hermitian operator of the spin current density 3 



31, 
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" V k (Ji 



(9) 



is a well-defined quantity (a tensor with nine compo- 
nents) only when the velocity operator is spin indepen- 
dent, as encountered in many metal spintronic devices. 45 
Such lack of physical justification for Eq. @ in SO cou- 
pled systems leads to an arbitrariness 1 ^ in the definition 



of the spin current density employed in recent intrinsic 
spin Hall studies, 46 thereby casting a doubt on the ex- 
perimental relevance of the quantitative prediction o 13 i 14 
for the spin Hall conductivity g s h = / E x computed 
as the linear response to the applied longitudinal elec- 
tric field E x penetrating an infinite SO coupled (perfect) 
semiconductor crystal. 

To obtain the spatial profiles of spin and charge cur- 
rent densities in finite-size samples of arbitrary shape at- 
tached to many probes, it is advantageous to represent 
the spin-dependent Hamiltonian and the corresponding 
charge and spin current density operators in the local or- 
bital basis ^SiAfiiilii 3 - For example, in such representation 
the Rashba Hamiltonian can be recast in the following 
form 9 



t aa ' ,r / / 
' ma ' mm' ^m'tr > 



(10) 



where hard wall boundary conditions account for confine- 
ment on the lattice L x x L y with the lattice spacing a. 
Here c mcr (c mcr ) is the creation (annihilation) operator of 
an electron at the site m = (ra x , m y ). 

While this Hamiltonian is of tight-binding type, its off- 
diagonal elements are non-trivial 2x2 Hermitian ma- 
trices t m / m = (t mm /)t in the spin space. The on-site 
potential e m describes any static local potential, such 
as the electrostatic potential due to the applied voltage 
or the disorder simulated via a uniform random variable 
e m G [—W/2,W/2]. The generalized nearest neighbor 
hopping t™ m , = (t mm f) aa f accounts for the Rashba cou- 
pling 



-t I a -itso&y (m = m' + e x ) 
-t Q I a + itso&x (m. = m! + e y ) ' 



(11) 



through the SO hopping parameter tso = a /2a (I s 
is the unit 2x2 matrix in the spin space). The di- 
rect correspondence between the continuous Eq. (J3J) and 
the lattice Hamiltonian Eq. (JTJJ) is established by using 
t — h 2 /2m*a 2 for the orbital hopping and by selecting 
the Fermi energy {Ep = — 3.8t in the rest of the pa- 
per) at which zero-temperature quantum transport takes 
place close to the bottom of the band Eb = — 4.0t o to 
ensure that injected quasiparticles have quadratic and 
isotropic energy-momentum dispersion which character- 
izes the Hamiltonians in effective mass approximation. 
Using the effective mass and conduction bandwidth of 
the semiconductor heterostructures measured in experi- 
ments, 47 one can interpret the parameters of the lattice 
Hamiltonian as having the typical values a ~ 3 nm for 
the lattice spacing while the Rashba SO hopping is of the 
order of tso/t Q ~ 0.01. 

The usage of the second quantized notation in Eq. El 
facilitates the introduction of Keldysh Green func- 
tion 42 expressions for the nonequilibrium expectation val- 
ues. 40,43 We imagine that at time t' = — oo the sample 
and the leads are not connected, while the left and the 
right longitudinal lead of a four-probe device are in their 
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own thermal equilibrium with the chemical potentials /i^ 
and hr, respectively, where hl — \ir + eV. The adia- 
batic switching of the hopping parameter connecting the 
leads and the sample generates time evolution of the den- 
sity matrix of the structured The physical quantities are 
obtained as the nonequilibrium statistical average (...) 
(with respect to the density matrix 42 at time t' = 0) 
of the corresponding quantum-mechanical operators ex- 
pressed in terms of c ma and c mcr . This will lead to the 
expressions of the type (cl na c ma t ), which define the lesser 
Green function 40,43 



?G.< 



m'm, a' a 



-f 

2m J_, 



dEG<, a (E). (12) 



Here we utilize the fact that the two-time correlation 



function [c ma (t) = e lHt / h c 



-iH 



G^^AU) = \ (tlAt')trnAt)) , (13) 

depends only on r = t — t' in stationary situations, so it 
can be Fourier transformed to energy 

G< m><T ,(r) = — h \ JEG< m ,^,(E)e iET / h , (14) 

which will be utilized for steady-state transport studied 
here. We use the notation where G^ m , is a 2 x 2 matrix 
in the spin space whose a a' element is G mm , acr , . 



A. Bond charge currents in SO coupled systems 

1. Bond charge- current operator 

Using the charge conservation expressed through the 
familiar continuity equation Eq. yields a uniquely de- 
termined bond charge current operator for quantum sys- 
tems described on a lattice by a tight-binding-type of 
Hamiltonian Eq. (|1U|). That is, the Heisenberg equation 
of motion 



dN n 
~dt 



1 

ih 



(15) 



for the electron number operator N m on site m, 



(16) 



leads to the charge continuity equation on the lattice 
dN„ 



dt 



(^m,m+e fc 

k=x,y 



i)=o. (17) 



This equation introduces the bond charge-current oper- 
ato r 34 i 39 J mm / which describes the particle current from 



site m to its nearest neighbor site m' (the 'bond' termi- 
nology is supported by a picture where current between 
two sites is represented by a bundle of flow lines bunched 
together along a line joining the two sites). 

Thus, the spin-dependent Hamiltonian Eq. (fTTTj) con- 
taining the non-trivial 2x2 matrix hopping defines the 
bond charge current operator J mm / = J m °m> which 

can be viewed as the sum of four different spin-resolved 
bond charge- current operators 



j 1 



e 

ih 



u m'a' L m'm 1 ma 



h.C. 



(18) 



where h.c. stands for the Hermitian conjugate of the 
first term. In particular, for the case of t™ m , being de- 
termined by the Rashba SO interaction Eq. ([TT]). we can 
decompose the bond charge current operator into two 
terms 



T / — T 

u mm' — rJ i 



kin 



fSO 



(19) 



which have transparent physical interpretation. The first 
term 



J, 



kin 



et Q 

ih 



At 



C w /^c mcr h.c. 



(20) 



can be denoted as "kinetic" since it originates only from 
the kinetic energy t Q and does not depend on the SO 
coupling energy tso, while the second term 



j 



so 



-—^ S mm> (m^m + e*) 

46^SO At / / i \ 

~^2-5 mm , (m =m + e 2/ ) 



4etso 

h 2 



((m'-m) x S mm /) (21) 



is the additional contribution to the intersite current flow 
due to non-zero Rashba SO hopping tso- Here we also 
introduce the "bond spin-density" operator 



cx(3 



(3Cmf3 



h.C. 



(22) 



defined for the "bond" connecting the sites m and m', 
which reduces to the usual definition of the local spin 
density operator for m' = m [see Eq. ([34)1 ]. 



2. Nonequilibrium bond 



current 



The formalism of bond charge currents yields the phys- 
ically measurable 38 local charge current within the sam- 
ple as the quantum- statistical average (...) (with respect 
to a density matrix that has evolved over sufficiently long 
time so that nonequilibrium state and all interactions are 
fully established) of the bond charge-current operator in 
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the nonequilibrium state^S^i^ 



^1 \^mm' 

-e r°° dE 

^ J-oo 2^ 



(23) 



m'm lT mm' .era 



,(E) 



t mm /G m , m f(T (E) ,(24) 



where we utilize Eq. (p"2|) to express the local charge cur- 
rent in terms of the nonequilibrium lesser Green function. 
The spin-resolved bond charge current in describes 
flow of charges which start as spin a electrons at the 
site m and end up as a spin a' electrons at the site m' 
where possible spin-flips a / a' (instantaneous or due 
to precession) are caused by spin-dependent interactions. 
The decomposition of the bond charge-current operator 
into kinetic and SO terms in Eq. (|19|) leads to a Green 
function expression for the corresponding nonequilibrium 

bond charge currents (|j mm /^ = ^mm'^ + ^mm') w ^ n 
kinetic and SO terms given by 



kin 



et, 



r /lfr Trs [ G ^(^)-G<, m (^)], 



(25) 



J 



so 



etso 



dE 



h J-oo 27 ™ 

xTr s [((m' - m) x a) z (G< m ,(£) + G<, m (i5))] . 

(26) 



Note, however, that kinetic term is also influenced by 
the SO coupling through G < . In the absence of the SO 
coupling, Eq. vanishes and the bond charge current 
reduces to the standard expression 42^^ The trace Tr s 
here is performed in the spin Hilbert space. Similarly, we 
can also obtain the nonequilibrium local charge density 
in terms of G < 



e y , { C ma C ma) 

e 
2iti 



/oo 

= — J^dETr s [G< m (E)], (27) 

which is the statistical average value of the corresponding 
operator Eq. (|16j). 

B. Bond Spin Currents in SO Coupled Systems 

1. Bond spin-current operator 

To mimic the plausible definition of the spin-current 
density operator in Eq. Q, we can introduce the 
bond spin-current operator for the spin-S'i component as 



the symmetrized product of the spin--| operator h&i/2 
(i = x, y, z) and the bond charge-current operator from 
Eq. (EJ) 

By inserting the hopping matrix t m / m Eq. of the 
lattice SO Hamiltonian into this expression we obtain its 
explicit expression for the Rashba SO coupled system 

^mm' = ~2i E {fm'(3 ( a i) (3 a ^ma ~ h.C. 



a(3 



tso 



iV mm ' (e* x (m' - m)) z , (29) 



which can be considered as the lattice version of Eq. 
Here N m m' is the bond electron-number operator 



^ = lj2 (i/ma + h.C.) , 



(30) 



which reduces to the standard electron-number operator 
Eq. Jini for m' = m. 



2. Nonequilibrium bond spin current 

Similarly to the case of the nonequilibrium bond charge 
current in Sec. II A 21 the nonequilibrium statistical aver- 
age of the bond spin-current operator Eq. (|2HJ) can be 
expressed using the lesser Green function G < as 



J, 



Ji 



■Si (kin) 



Si (kin) 



j; 



^ (so) 



(31) 



00 dE 

oo 27r 



xTr[a, (G< m ,(£)-G<, m (£))],(32) 

fSi(so)\ / / / \ \ ^so f°° dE 

xTr[G< m ,(E) + G< >m (E)]. (33) 

Here we also find two terms which can be interpreted 
as the kinetic and the SO contribution to the bond spin 
current crossing from site m to site m'. However, we 
emphasize that such SO contribution to the spin-S^ bond 
current is identically equal to zero, which simplifies the 
expression for this component to Eq. (|52J) studied in the 
rest of the paper as the primary spin current response in 
the spin Hall effect. 



3. Local Spin density and its continuity equation 

The local spin density in the lattice models is deter- 
mined by the local spin operator S m = (S^, ^m, *§m) a ^ 
site m defined by 



h 



(34) 



a(3 



6 



The Heisenberg equation of motion for each component 
Si (i = y, z) of the spin density operator 



dt 



1 

ih 



S l m ,H 



(35) 



can be written in the following form 



dt 



k=x,y 



fS t 



^mS (36) 



where J^m' ^ s tne bond spin-current operator given by 
Eq. (|28j) so that the second term on the lefthand side 
of Eq. (|36|) corresponds to the "divergence" of the bond 
spin current on site m. Here, in analogy with Eq. 0, we 
also find the lattice version of the spin source operator 
F^* whose explicit expression is given by 



m 


_ *SO / f Sz 

t V m ' m+e ^ 


+ J s * \ 

' u m—e x ,m J J 


(37) 




*SO / f Sz 
^ \^m,m+e y 


+ J 8 ' \ 

1 ' m—e y .m J i 


(38) 


F Sz 


_ ^so / f Sx 


fs x 

' me, ,m 





, J S y 

1 ' m.m+e ( . 



,J m e,. 



• (39) 



The non-zero term F*j* on the righthand side of the 
spin continuity equation Eq. is a formal expression, 
within the framework of bond spin current formalism, of 
the fact that spin is not conserved in SO coupled sys- 
tems where it is forced into precession by the effective 
momentum-dependent magnetic field of the SO coupling. 
The fact that the bond spin current operator Eq. (|28|) ap- 
pears in the spin continuity Eq. as its divergence im- 
plies that its definition in Eq. (|28j) is plausible. However, 
the presence of the spin source operator F*j* reminds us 
that such definition cannot be made unique^ unlike the 
case of the bond charge current which is uniquely deter- 
mined by the charge continuity Eq. ([T7j). 

If we evaluate the statistical average of Eq. (|35J) in a 
steady state (which can be either equilibrium or nonequi- 
librium), we obtain the identity 



' m.m+e/ / \ u m e /, . m 



)={F%). (40) 



k=x,y 

In particular, for the spin-5^ component we get 

E 



k=x,y 



■ j>, \ _ / jS z \\ 

vr 1 vn-Uo, / \ fJ e / . m / J 



tso 

to 



.,m+e fc 



fSk 



k=x,y 



),(41) 



which relates the divergence of the spin-S^ current (left- 
hand side) to the spin-source (righthand side) determined 
by the sum of the longitudinal component of the spin- 5^ 



current and the transverse component of the spin-5^ cur- 
rent. 

Since no experiment has been proposed to measure lo- 
cal spin current density within the SO coupled sample, 
defined through Eq. or its lattice equivalent Eq. (|31[) , 
we can obtain additional well-defined and measurable in- 
formation about the spin fluxes within the sample by 
computing the local spin density 

H 



/CO 
dE E ^G< m>a0 (E) 

-°° a,/3=T,J. 

= — J JETv s {aG< m (E)}. (42) 

Plotting of this quantity for Rashba SO coupled systems 
will provides us with insight into the stationary spin flow 
in the nonequilibrium steady transport state. 



C. Spin-resolved Landauer-Keldysh Green 
functions for finite-size mesoscopic devices 

The formalism discussed thus far does not depend on 
the details of the external driving force which pushes the 
system into a nonequilibrium state. That is, the system 
can be driven by either the homogeneous electric field 
applied to an infinite homogeneous 2DEG or the volt- 
age (i.e., electrochemical potential) difference between 
the electrodes attached to a finite-size mesoscopic sam- 
ple. For example, in the latter case, the external bias 
voltage only shifts the relative chemical potentials of the 
reservoirs into which the longitudinal leads (employed to 
simplify the boundary conditions) eventually terminate, 
so that the electrons do not feel any electric field in the 
course of ballistic propagation through clean 2DEG cen- 
tral region. The information about these different situa- 
tions is encoded in the lesser Green function G < . 

Here we focus on experimentally accessible spin Hall 
bridges^ where finite-size central region (C), defined on 
the LxL lattice, is attached to four external semi-infinite 
leads of the same width L. The leads at infinity termi- 
nate into the reservoirs where electrons are brought into 
thermal equilibrium, characterized by the Fermi-Dirac 
distribution function f(E — eV p ), to ensure the steady- 
state transport — in such Landauer setup^i current is lim- 
ited by quantum transmission through a potential profile 
while power is dissipated non-locally in the reservoirs. 
The voltage in each lead of the four-terminal spin Hall 
bridge is V p (p = 1, . . . , 4) so that the on-site potential 
e m within the leads has to be shifted by eV p . 

The spin-dependent lesser Green function G < defined 
in Eq. (|13[) is evaluated within the finite-size sample re- 
gion as a 2L 2 x 21? matrix in the site(g)spin space through 
the spin-resolved Keldysh equation for matrices 9 



G<(E) = G(E)J: < (E)G\E), 



(43) 



7 



which is valid in this form for steady-state transport when 
transients have died away™ Within the effective single- 
particle picture, the retarded Green function can be ob- 
tained by inverting the Hamiltonian 



G(E) 



EIc -Hc-U m -J2 s p(£ " eV p) 



where the self-energies 



(44) 



£<(£) = iJ2r p (E-eV p )f(E-eV p ), 



[V P (E) - St (E)} , 



i0+) I p - H 



lead! 



(45) 

(46) 
H pC , (47) 



are exactly computable in the non-interacting electron 
approximation and without any inelastic processes tak- 
ing place within the sample. They take into account the 
"interaction" of the SO coupled sample with the attached 
leads in the Landauer transport setup where they gen- 
erate finite that electrons spends within the 2DEG be- 
fore escaping through the leads toward the macroscopic 
thermalizing reservoirs. Here Ic is the 21? x 21? iden- 
tity matrix, while I p is the identity matrix in the infinite 
site(g)spin space of the pth lead, and we use the following 
Hamiltonian matrices 



(■^■c) mm' ,aa' 

(H' ead ) , , 
(H P c) r 



/ mm' ,aa' 



(lnurl^lW), (m,m'eC), 
(l mcr | H |l m v) , (m,m'ep), 
{l m a\H\l m , a ,) , (mGp.m'eC), 

(48) 



with |l mcr ) being a vector in the Fock space (meaning 
that the occupation number is one for the single particle 
state \m<j) and zero otherwise) and H is the Hamiltonian 
given in Eq. (fTHjl . 

In the general case of arbitrary applied bias voltage, 
the gauge invariance of measurable quantities (such as 
the current-voltage characteristic) with respect to the 
shift of electric potential everywhere by a constant V, 
V p — > V p + V and U m — > U m +V, is satisfied on the proviso 
that the retarded self-energies J} p (E — eV p ) introduced 
by each lead depend explicitly on the applied voltages at 
the sample boundary, while the computation of the re- 
tarded Green function G(E) has to include the electric 
potential landscape U m within the sample 48 [which can 
be obtained from the Poisson equation with charge den- 
sity Eq. (|77j)]. However, when the applied bias is low, 
so that linear response zero-temperature quantum trans- 
port takes place through the sample [as determined by 
G(Ef)], the exact profile of the internal potential be- 
comes irrelevant. 34,49 



II. SPATIAL DISTRIBUTION OF LOCAL SPIN 
CURRENTS AND SPIN DENSITIES IN 
BALLISTIC FOUR-TERMINAL RASHBA SO 
COUPLED NANOSTRUCTURES 

Under the time reversal transformation t — > — t, the 
mass, charge, and energy are even while the velocity and 
Pauli matrices are odd. Thus, the charge currents are t- 
odd and must vanish in thermodynamic equilibrium (in 
the absence of magnetic field). On the other hand, the 
spin current-density operator Eq. (0), which is the prod- 
uct of two t-even quantities, can have its expectation val- 
ues to be non-zero even in thermodynamic equilibrium. 
This has been explicitly demonstrated 1 ^ for its expecta- 
tion values in the eigenstates of an infinite clean Rashba 
spin-split 2DEGi2 

To investigate possible patterns of such equilibrium lo- 
cal spin currents in mesoscopic finite-size devices we plot 
in Fig. [Ha) the spatial distribution of the bond spin cur- 
rents, carried by the whole Fermi sea, in a four-terminal 
ballistic device with no impurities where all leads are kept 
at the same potential. Although we find non-zero local 
spin currents, they do not transport any spin since the 
total spin current, obtained by summing the bond spin 
currents over an arbitrary transverse cross section of the 
device 



^trans( m 2/) y(m x ,m y )(m x ,rn y + l)/ ' ( 49 ) 

or over any longitudinal cross section 



(m x ,m y )(m x + l,m y ) 



(50) 



is identically equal to zero It rans (m y ) = ^iong( m ^) = 0- 
This picture also provides a direct microscopic proof 
that no equilibrium total spin currents, conjectured in 
Ref. |33, can actually appear in the leads of an unbiased 
V p = const, mesoscopic device in thermodynamic equilib- 
rium, as demonstrated recentlyi2i£S within the Landauer- 
Biiittker approach which operates only with the total 
spin and charge currents in the leads. 

In Fig.. [IJb) we apply low (i.e., linear response, see 
Sec. [TVJ) bias voltage eV = 10'% <C (E F - E b ) = 0.2t o 
between the longitudinal leads and integrate expression 
Eq. (|32|) from the bottom of the band to the chemical 
potential Ep + eV/2 of the left reservoir. In contrast 
to the equilibrium spin current density from Fig. QJa), 
the vortex pattern is now distorted and non-zero total 
spin current It rans (m y ) ^ in Eq. emerges in the 
transverse direction, as expected in the phenomenology 
of the spin Hall effect. 

One of the highly unconventional features of the in- 
trinsic spin Hall current is its dependence on the whole 
SO coupled Fermi sea^iii even when infinite system is 
driven out of equilibrium by the applied external electric 
field thereby limiting the charge transport and the ex- 
trinsic spin Hall response to the Fermi level through the 
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FIG. 1: (Color online) The spatial distribution of the local spin current in ballistic four-terminal bridges where the central 
2DEG region, with the Rashba SO coupling tso = 0.1t o and the corresponding spin precession length Lso = 7rt a/2tso ~ 15.7a 
(typically a c± 3 nm), is attached to four ideal (tso = 0) leads. The magnitude of the bond spin current is proportional to the 
length of the arrow. The device is in equilibrium in (a), and out of equilibrium in (b), (c), and (d) due to the applied bias 
voltage eV = 10 _3 t o which drives the linear response longitudinal charge current and the transverse spin Hall current induced 
by its passage through the SO coupled region. The local spin current in (b), which is "carried" by all states from — 4t (band 

bottom) to Ef + eV/2 (Ef = — 3.8t ), is the sum of the equilibrium (persistent) spin current ^mm^} m ( c )' carr i e d by the 

fully occupied states from — 4t Q to Ef — eV/2, and the nonequilibrium (transport) spin current (j^^t^^ in (d) carried by 

the partially occupied states around the Fermi energy from fiR — Ef — eV/2 (electrochemical potential of the right reservoir) 
to hl = Ef + eV/2 (electrochemical potential of the left reservoir). 



nonequilibrium part of the distribution function. 21 How- 
ever, such property appears to be alien to the spirit of 
Fermi liquid theory where transport involves only quasi- 
particles with energies within /c#T of the Fermi level. 



While the spin-currents crossing the transverse bonds 
in Fig. IHb) are apparently carried by the whole Fermi 
sea, we now separate the integration in Eq. (|32|) for S z 
bond spin current into two parts 



E F -eV/2 E F +eV/2 



Jt m ) = | / ^Tr [*, (G< m ,(E) - G<, m (E))] + | J ^Tr [&, (G< m ,(£) - G<, m (E))} 

E h E F -eV/2 

&) + (&)■ (51) 
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The states from the band bottom E\> to Ep — eV/2 are 
fully occupied, while states in the energy interval from 
the electrochemical potential Ep — eV/2 (e < 0) of the 
right reservoir to the electrochemical potential Ep + eV/2 
of the left reservoir are partially occupied because of the 
competition between the left reservoir which tries to fill 
them and the right reservoir which tries to deplete them. 

The profile of the first term ^ j^,^ ^ in Eq. (|51j) is shown 

in Fig.[IJc), while the spatial profile of the second term, 

representing the local spin current (^J^m'^^ carried by 

the states around the Fermi energy is shown in Fig.^d). 

The spatial distribution of the microscopic spin cur- 
rents in Fig. EJc) is akin to the vortex-like pattern of 
bond spin currents within the device in equilibrium in 
Fig. [Ha), and, therefore, does not transport any spin be- 
tween two points in real space. Thus, Fig. convincingly 
demonstrates that the non-zero spin Hall flux through 
the transverse cross sections in Fig. QJb) and Fig. QJd) 
is due to only the wave functions (or Green functions) 
at the Fermi energy (as T — > 0), in accord with the gen- 
eral paradigms of the Landau's Fermi liquid theory where 
transport quantities are expected to be expressed as the 
Fermi-surface property^ii 

We recall here that similar situation appears in charge 
transport in an external magnetic field where equilibrium 
(or persistent) current density^ or bond charge currents 
in the lattice formalism f^LAi is non-zero even in unbi- 
ased devices (all leads at the same potential) in thermal 
equilibrium due to breaking of time-reversal invariance 
by the external magnetic field. However, such circulat- 
ing or diamagnetic currents carried by the Fermi sea, 
which in Landauer-Keldysh formalism can be subtracted 
by separating the integration^! in a fashion similar to 
our Eq. (|5l|) . do not contribute to the net charge trans- 
port (i.e., to the total charge current measured in ex- 
periments) through any cross section of the device^^i 
Thus, early "Fermi sea" expressions for the linear re- 
sponse transport coefficients in, e.g., the quantum Hall 
effect theory 34 or in the anomalous Hall effect theory^ 51 
were eventually recast in terms of the Fermi surface de- 
termined quantities. Similarly, Fig. E] demonstrates that 
spin Hall current carried by the "bulk" of the Fermi sea, 
which is equilibrium and does not really transport spin 
between two points in space, must be subtracted^ within 
any sensible theory for the intrinsic spin Hall conductiv- 
ity defined in the thermodynamic limit of macroscopic 
systems. 

One of the basic tests for theories of the spin Hall effect 
is to predict the direction (i.e., the sigrAALffi) of the spin 
Hall current or the corresponding nonequilibrium spin 
Hall accumulation deposited by such current on the lat- 
eral boundaries of experimental devices i£ The nonequi- 
librium (linear response) spin Hall current in Fig. H^d) 
flows from the top to the bottom transverse lead be- 
cause the spin- 1 electrons are deflected to the right. This 
feature can be understood using the semiclassical pic- 
ture based on the SO force operator generated by the 



Rashba Hamiltonian Eq. Q of the finite-size 2DEG 

- 2a 2 m* //s N ^ dV con Ay) , . 

Fso = — ^— (p x z) <g> a* ^gp^Y* ( 52 ) 

and its expectation value in the spin-polarized wave 
packet state <g> ||). This apparently simple picture 
also explains why the transverse spin Hall current density 
bends toward the right in Fig. while passing through 
the SO coupled region. 

However, this expectation value (i.e., the SO "force") 
oscillates along the sample due to the precession of the 
deflected spin in the effective Rashba magnetic field 
which is nearly parallel to the ?/-axis because of the trans- 
verse confinement effects. 52 In ballistic strongly coupled 
SO structures such a 2 -dependent SO "force", which os- 
cillates on the mesoscale set by the spin precession length 
Lso = 7rh 2 /2m*a = 7rt a/2tso (on which the spin pre- 
cesses by an angle 7r, i.e., the state | j) evolves into ||)), 
will lead to a change of the signi 5 - of spin Hall current 
as a function of the system size L/L^o- Also, since the 
mesoscopic spin Hall effect sensitively depends on the 
measurement geometry^ 5 the sign of the transverse spin 
Hall current can change when non-ideal leadsii are at- 
tached to the sample. 

To highlight mesoscopic features^ (such as the effect 
of the measuring geometry and the properties of the at- 
tached probes) of the spin Hall effect in multiterminal 
ballistic SO coupled structure^ we plot in Fig. [21 the 
spatial profile of microscopic spin currents for the four- 
terminal bridge attached to the leads where Rashba SO 
coupling is switched on adiabatically (via a linear func- 
tion) within the finite region of the leads adjacent to the 
2DEG sample in the center of the device. In this measur- 
ing setup, the reflection 5 ^ at the interface separating zero 
and non-zero SO coupling regions is greatly suppressed, 
thereby enhancing the spin Hall current [as encoded by 
longer arrows in the profiles of Fig. [21 when compared to 
Fig. Del)]. 

Since local spin current within the central region is 
not conserved, as manifested by the total spin current 
changing magnitude between different transverse cross 
sections (see Sec. II VI for total spin current profiles) sep- 
arated by long > Lso> the s P a tial profiles of the local 
spin current density appear not to be directly measur- 
able. 3 However, the nonequilibrium spin density (i.e., the 
related spin magnetization) is a well-defined and measur- 
able quantity. 1,2 Therefore, we plot the spatial distribu- 
tion of the stationary flowing l^S^j Eq. spin density 

in Fig. [2 to gain an additional insight into the dynamics 
of the spin Hall quantum transport. These microscopic 
picture convincingly demonstrate how spin-j and spin-j 
densities flow in opposite transverse directions through 
the attached leads, without any net charge flow in the 
transverse direction, thereby giving rise to a pure trans- 
verse spin Hall current. When we reverse the direction 
of the longitudinal charge current (by reversing the bias 
voltage V — > —V), the transverse spin current and spin 
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FIG. 2: (Color online) (a) The spatial distribution of the nonequilibrium local spin current (^J^m'^^) ^ n ballistic four-terminal 
spin Hall bridges where the central finite-size 2DEG (of size 8a x 8a) with the Rashba SO coupling tso = 0.1t o (Lso ~ 15.7a) 
is attached to four leads containing a region (adjacent to the sample) of length 8a within which the SO coupling is switched on 
adiabatically (using a linear function) from tso = to tso = 0.1t o . Panel (b) shows the spatial profiles of the stationary flow 

of the local spin density (^S^, while panel (c) plots the spatial distribution of the local charge current density (^J mm >^. Upon 

reversing the bias voltage V — > —V driving the linear response (eV = 10 _3 t o ), spatial profiles in the left column turn into the 
profiles of the right column where both the total longitudinal charge current and the total transverse spin Hall current change 
their direction?- 
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FIG. 3: (Color online) (a) The spatial distribution of the nonequilibrium local spin current (j^m'^ / an< ^ O 3 ) ^ ne flowing 

spin density i n ballistic bridges consisting of four ideal leads (tso = 0) attached to a finite-size 2DEG with the Rashba 

SO coupling tso = 0.1t o and the corresponding spin precession length Lso ~ 15.7a. The 2DEG central region is of the size 
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densities flip their sign, as exploited in experiments to 
confirm strong signatures of the spin Hall effect. 2 

Thus, in contrast to the arguments 21 advocating im- 
possibility of spin transport and accumulation by mech- 
anisms driven solely by the intrinsic SO coupling terms 
in the effective Hamiltonian of spin-split semiconductors, 
Fig. demonstrates that spin Hall effect originating in 
ballistic multi-terminal devices (without necessity for im- 
purity induced effects) generates genuine nonequilibrium 
spin flux that can be used for spin injection and spin- 



tronics applications A 

Another mesoscale driven property of the spin Hall 
effect in ballistic nanostructures is its finite-size scal- 
ing being governed by the processes on the spin preces- 
sion length Lso' One can differentiate the "mesoscopic" 
regime L < Lso where spin Hall current oscillates with 
increasing 2DEG size (changing sign with increasing size 
L x L of the 2DEG for strong enough Rashba SO coupling 
^so ^ 0.04£ o ); and the "macroscopic" regime L ^> Lso 
where it saturates at some average value. 15 While the 
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FIG. 4: (Color online) (a) The spatial distribution of the nonequilibrium local spin current (j^^ eq ^) an d (b) the flowing spin 

density (^S^ in the disordered 2DEG of the size L — 60a with static spin-independent impurities which set the mean free path 

£ = 7a. The Rashba SO coupling within the 2DEG region is tso = 0.1£ o (tso = in ideal the leads) and the corresponding 
D'yakonov-Perel' spin relaxation length 5 is Lso ~ 15.7a. The applied bias voltage eV = 10 _3 t o between the longitudinal leads 
drives the linear response longitudinal charge transport in the diffusive regime (L ^> £). Panels in the left column correspond 
to a single sample with specific configuration of impurities, while panels in the right column show disorder- averaged profiles 
over an ensemble of 100 different impurity configurations. 



spatial profiles of local spin currents and spin densities 
are easy to interpret for L < Lso ( as in Fig. Fig. [31 
divulges how they become increasingly more intricate in 
the samples of the size Lso x ^so (f° r which the spin 
Hall current reaches its maximum 15 ), or in the macro- 
scopic regime L ^> Lso- Nevertheless, the total spin 
Hall current exists only in the transverse leads (the sum 
of bond 5^-spin currents over the cross section of longitu- 
dinal leads is equal to zero) with its magnitude and sign 
being identical to the terminal spin currents obtained 



within the Landauer-Biittiker formalism 15 (see Sec. II VI 
for quantitative comparison). 



III. BULK VS. EDGE LOCAL SPIN CURRENTS 
IN DISORDERED FOUR-TERMINAL RASHBA 
COUPLED NANOSTRUCTURES 

A surprising feature of the early predictions for the 
intrinsic spin Hall effect in hole-doped 13 or electron- 



13 



doped 14 infinite homogeneous SO coupled semiconductor 
systems is apparent insensitivity of the spin Hall conduc- 
tivity <j s h on the mean free path and relaxation rates. 21 
This has provoked intense theoretical scrutiny of the ef- 
fects which spin-independent scattering off static impuri- 
ties 29 imposes on the intrinsic spin Hall current, leading 
to the conclusion that, in fact, cf s h — ► vanishes 23-25 for 
arbitrary small disorder in any model with SO coupling 
linear in momentum (such as the Rashba Hamiltonian of 
the spin-slit 2DEG) due to accidental cancellations. In 
the general case, where SO coupling contains higher or- 
der momentum terms jSiffi cr sH can be resilient to sizable 
disorder strengths. 54 Moreover, an influential conjecture 
has emerged from the studies of the disorder effects on the 
intrinsic spin Hall current — macroscopic inhomogeneities 
facilitate spin currents^ 6 so that transverse edge spin 
current J* could emerge near the sample-electrode in- 
terface, 27 even in systems where it is expected to be de- 
stroyed in the bulk^ 3 - 

However, the quantitative support for the picture of 
edge spin Hall currents is based on the analysis of semi- 
classical diffusive transport through a rather abstract 
structure, where 2DEG infinite in the transverse direc- 
tion is attached to two massive electrodes in the longi- 
tudinal direction, 27 that does not correspond to any ex- 
perimentally realizable device. On the other hand, the 
presence of SO couplings makes the dynamics of trans- 
ported spin in experimentally relevant confined struc- 
tures strongly dependent on the properties of their inter- 
faces, boundaries, and the attached electrodes, 53 - even for 
semiclassical spatial propagation of charges which carry 
spins evolving according to quantum dynamical laws^& 
For example, heuristic arguments based on the Keldysh 
formalism applied to an infinite two-terminal structure 
(lacking actual lateral edges) of Ref. |2j] suggest that 
nonequilibrium spin Hall accumulation (S^) ^ will ap- 
pear only in the four corners at the lead/2DEG interfaces 
(due to Jy 7^ existing within a spin relaxation length 
Lso wide region around such interfaces), in contrast to 
the Keldysh formalism applied to finite-size 2DEG in the 
Landauer two-terminal setup where non-zero spin accu- 
mulation (with opposite sign on the two lateral edgesiS) 
is found along the whole lateral edge. 9,10 

Within the formalism of bond spin currents of Sec. [I] 
(which represent the lattice version of the same quan- 
tity Jy studied in Ref. I27T) these issues can be resolved 
through the exact evaluation of the retarded and lesser 
Green functions a non-interacting particle propagating 
through the random potential in finite-size multiterminal 
SO coupled device. Thus, we plot in Fig. 0| the spatial 
distribution of the local spin currents and spin densities 
for a single disordered Rashba SO coupled 2DEG, as well 
as their disorder- averages (which are to be compared with 
the analysis based on the diffusion equation 27 ). The dis- 
order strength is tuned to ensure the diffusive transport 
regime £ <C L, while the magnitude of the spin Hall cur- 
rent in the leads at this concentration of impurities is still 
about 80% of its maximum value set in the clean limit 
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FIG. 5: (Color online) The longitudinal profile of the disorder- 
averaged (over 100 samples) nonequilibrium bond spin cur- 
rents, whose full two-dimensional profile is shown in the right 
panel of Fig. 21 a), across the transverse cross sections in the 
ideal bottom lead (m y — 5) and inside the diffusive 2DEG 
(m y = 20 and m y = 40) with static spin-independent im- 
purities. Note that the length of the m y — 5 cross section 
within the bottom transverse lead is 60a, while the length of 
m y — 20 and m y — 40 transverse cross sections is 80a (60a 
inside the 2DEG + 10a within each of the two longitudinal 
leads) . 



for the same four-terminal nanostructure^ 5 - We do not 
find any evidence for the confinement of spin fluxes near 
the 2DEG/longitudinal- leads interfaces. The conclusion 
based on the spatial profile of microscopic spin currents 
[i.e., the "distribution of arrows" in Fig. Ufa)] is further 
corroborated in Fig. by the one-dimensional longitudi- 
nal profiles of the bond spin currents on different cross 
sections within the SO coupled sample and in the leads. 
Thus, both Fig.^land Fig. [5| suggest that precessing spins 
will propagate through the bulk of the diffusive Rashba 
SO coupled 2DEG in semiconductor heterostructure. 



IV. TOTAL SPIN HALL CURRENT VS. LOCAL 
SPIN HALL CURRENT: LANDAUER-BUTTIKER 
VS. LANDAUER-KELDYSH PICTURE 

The traditional charge transport experiments measure 
total current I and the conductance I = GV relating 
the charge current to the voltage drop V, rather than 
local current density j and the conductivity j = aE re- 
lating it to the externally applied electric field (note also 
that in ballistic transport or quantum-coherent transport 
through a diffusive conductor conductivity ceases to exist 
as a local quantity 34 ). Since realization of the total pure 
spin currents has been detected experimentally in op- 
tical pump-probe experiments, 57 and several theoretical 
schemes are proposed to detect them indirectly via var- 
ious electrical measurements jLSLS we focus in this sec- 
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FIG. 6: (Color online) The total spin current Itvansi^y) on 
transverse cross sections across the ideal bottom and top 
transverse leads (with no SO coupling) and the Rashba SO 
coupled 2DEG with tso = 0.1t o (Lso « 15.7a). The linear 
response longitudinal charge transport (at zero temperature) 
is driven by the applied bias voltage eV = 10 _3 t o between the 
longitudinal leads of the four-terminal structures where un- 
polarized electrons are injected from the left longitudinal lead 
at the Fermi energy Ef = — 3.8t . The total pure transverse 
spin current in panel (a) is obtained by summing the bond 
spin currents in Fig. dd) over a cross section m y ; similarly, 
panel (b) and (c) corresponds to such sums for the profiles in 
the left and right columns of Fig. |3] respectively. 



tion on the properties of total spin Hall current /trans( m 2/) 
on different transverse cross sections of four-terminal de- 
vices, which is obtained from Eq. ((49)) by summing the 
bond spin currents. 

As shown in Fig. the total pure spin current in the 
transverse leads, obtained by summing the nonequilib- 



FIG. 7: (Color online) The longitudinal profile of the nonequi- 
librium bond spin currents along the transverse cross sections 
m y in the ideal bottom lead (solid curves) and inside the 
SO coupled ballistic 2DEG (dashed and dotted curves). The 
Rashba SO coupling is tso = 0.1t o (Lso ~ 15.7a) within 
the central 2DEG region and tso = in both the longitudi- 
nal and the transverse leads of the four-terminal bridge. The 
sum of the longitudinal profiles gives the total transverse spin 
Hall current in Fig. |S| Note that profile in panel (a) corre- 
sponds to bond spin current magnitude represented by arrows 
in Fig. Q^d) and, similarly, panels (b) and (c) corresponds to 
the left and right columns in Fig. |3| respectively. 



rium bond spin currents (^J^m^^ over an arbitrary 
transverse cross section of the ideal leads (where SO cou- 
pling vanishes), flows through them in a conserved fash- 
ion, i t s rans (m y ) = const, for any m y G lead. However, the 
same summation over the transverse cross sections within 
the 2DEG yields a quantity which is not conserved, ex- 
cept on the short length scales <C £so This is due to the 
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We demonstrate in Fig. |HI that the spin Hall conduc- 
tance obtained from the Landauer-Biittiker formalism 
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FIG. 8: (Color online) The spin Hall conductance of a 
Rashba SO coupled four-terminal 2DEG, obtained from the 
Landauer-Biittiker linear response transmission formalism as 
G z }^^ — limy^o 1%/V [where If denotes the terminal spin 



current in the top transverse lead in Fig. OH a) 



Hall conductance G 



*(LK) 



the spin 



A S rans( m y) /V obtained within the 



Landauer-Keldysh formalism by summing the bond spin cur- 
rents over the transverse cross section in the top lead of struc- 
ture in Fig. Eta) for low bias voltage eV = 10 _3 t o or high bias 
voltage eV = 0.1t o applied between the longitudinal leads. 
The size of the central 2DEG region is 15a x 15a « Lso x Lso, 
for which the spin Hall conductance G z }^^ (L, tso) reaches 
maximum when increasing the sample size at fixed SO cou- 
pling^ 



fact that, e.g., injected eigenstate | j) of & z will precess 
in the effective magnetic field of the Rashba SO cou- 
pling (along the y-axis), thereby changing the amplitude 
of the spin current measured with respect to the z-axis 
as the spin quantization axis. Additional quantitative 
information about the microscopic details of spin fluxes 
is provided by Fig. [7| which plots the one-dimensional 
longitudinal profiles of the bond spin currents over the 
selected transverse cross sections (in the leads and in the 
2DEG sample) cutting through the full spatial distribu- 
tions of Fig. [l] and Fig. 03 Note that the sum of these 
longitudinal profiles yields the corresponding total spin 
current at the cross section m y in Fig. El 

The total spin currents in the leads in the linear re- 
sponse regime V — > can also be calculated using the 
spin-dependent Landauer-Biittiker scattering formalism 
for spin currents in multiprobe geometries iU&Si In this 
formalism, one has to compute the spin-resolved trans- 
mission matrices connecting the spin-polarized asymp- 
totic scattering states in semi-infinite ideal leads at- 
tached to the sample, 53 which requires the knowledge 
of only the elements of the retarded Green function ma- 
trix Eq. (|44|) between the sites of the sample which are 
adjacent to the leads (similar expressions in terms of 
only the retarded Green function between the sample 
edges are obtained for the total charge currents^^ 3 -) . 



z(LB) 



limy ^qI^/V (where 2 denotes the top trans- 



G 

verse leacU 5 -) is almost identical to the spin Hall con- 
ductance G z s ff K ^ = ^trans( m ?/)/^ obtained by summing 
the bond spin currents over a cross section m y in the top 
transverse ideal lead, on the proviso that the applied bias 
voltage is small eV <C Ep and current is carried only by 
the states at the Fermi level. This result further justi- 
fies the introduction of nonequilibrium bond spin current 
formula Eq. (|H2[) for the z-component of spin. 



V. CONCLUDING REMARKS 

In conclusion, we have demonstrated how to define the 
bond spin current, describing the spin flux across a sin- 
gle bond between two sites of the lattice model of an SO 
coupled semiconductors, and evaluate it in terms of the 
Keldysh nonequilibrium Green functions for to the Lan- 
dauer setup where finite-size sample is attached to many 
semi-infinite ideal leads to form a theoretical model of ex- 
perimentally accessible spin Hall bridges. Although spin 
current is not conserved within the SO coupled region 
(i.e., on the length scales comparable to the spin preces- 
sion length Lso on which the SO coupling manifests it- 
self), the microscopic spin fluxes are nearly conserved on 
short scales. Thus, the bond spin currents make it pos- 
sible to obtain their spatial distribution by following the 
dynamics of transported spin on the scale of the lattice 
spacing a <C Lso- Such profiles of (the lattice version of) 
local spin current density, together with stationary flow 
profiles of physically transparent (and measurable) local 
spin densities, allows us to demonstrate microscopic de- 
tails of how pure transverse spin Hall current emerges in 
clean Rashba SO coupled 2DEG through which the un- 
polarized longitudinal charge current flows ballistically 
(where electrons do not scatter off impurities and do not 
feel the electric field). These spatial profiles are highly 
dependent on whether the 2DEG size is smaller or greater 
than the spin precession length, and can be affected by 
non-trivial measurement geometries, as discussed in the 
theory of the mesoscopic spin Hall effect. 15 The bond spin 
current within the bulk of the 2DEG is also resilient to 
weak disorder so that spin fluxes remain non-zero in the 
bulk of the sample and are not localized near the edges 
of the diffusive Rashba spin-split 2DEG. 

Using the bond spin current we explicitly demonstrate 
that nonequilibrium total spin current can be carried only 
by the states around the Fermi energy, while the Fermi 
sea contributes to local persistent spin currents which, 
however, do not transport any spin through a given cross 
section. The nanometer scale details of the spin Hall flow 
in multiterminal mesoscopic (quantum-coherent) struc- 
tures convincingly show that the intrinsic to the crystal 
SO couplings can be used to generate spin fluxes, spin 
accumulation, and ultimately be employed to construct 
all-electrical spin injectors which do no require any fer- 
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